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Collagen fibrils, cable-like assemblies of long biological molecules, the so-called triple 
helices, are dominant components of connective tissues. Their determinant morphological 
and functional roles motivated a large number of studies concerning their formation and 
structure. However, these two points are still open questions and, particularly, that of 
the lateral assembly of the triple helices which is certainly dense but not strictly that 
of a well-ordered molecular crystal. We examine here the geometrical template provided 
by the algorithm of phyllotaxis which gives to each element of an assembly of points or 
parallel rods the most homogeneous and isotropic dense environment in a situation of 
cylindrical symmetry. The scattered intensity obtained from a phyllotactic distribution of 
triple helices in collagen fibrils presents features which could contribute to the scattering 
observed along the equatorial direction of their X ray patterns. Following this approach, 
the aggregation of triple helices in fibrils should be considered within the frame of soft 
condensed matter studies rather than that of molecular crystal studies. 
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1. Introduction 

Collagen fibrils are the major constituent of connective tissues such as bone, carti- 
lage, myofibrils, ligament, skin, cornea. Typical fibrils of type I collagen are shown 
in figure [TJ Their properties and the rich polymorphism of their associations under 
different physico-chemical conditions allow the building of tissues adapted to the 
very diverse roles expected from them. Unfortunately, they can also underly the 
development of crippling deformities associated with rheumatic diseases and con- 
genital defects in the tissues. The structure of collagen fibrils has therefore been the 
subject of many studies but is not yet totally solved. It is now commonly accepted 
that the formation of type I collagen fibrils proceeds along two steps as shown in 
figure [2] left-handed simple helices of collagen molecules build right handed triple 
helices which in turn build the fibrilfl If the first step is now well described and 

^An intermediate step in which five triple helices would associate to build microfibrils was proposed 
earlier but is now questioned. 
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Fig. 1. Some fibrils of type I collagen observed with electron microscopy. Their diameter is of 
about 150 nm, the period of their longitudinal striations is of 67 nm (from J. Gross). These 
morphological features are rather constant whatever the conditions, in vivo as well as in vitro. 



understoodij, the second, the lateral organization of triple helices, is still subject of 
studies which can be summarized as follows. 

The main information concerning the lateral organization of triples helices in 
fibrils is provided by X-ray scattering spectra whose intensities are directly related 
to the square of the Fourier transform of the organization. Typical X ray patterns are 
shown in figure [3] one from HulmeJ^, a drawing showing the results of backgro und 
eliminations to put in light weak localized features as proposed by Wess and OrgePEl 
and another pattern from DouceM 

These examples show that all the intensity is concentrated in the wave vector 
range k < 5nm^^ only. In all cases, an intense scattering is observed for k w 5nm^^ 
which is to be related to the packing distance between triple helices as, if they 
were considered as straight parallel rods with a 1 nm diameter packed at a distance 
ff 1.3 nm along a hexagonal lattice, this should be the position of the first Bragg 
peak. In all cases also an important fan-like background scatter, including equatorial 
central scatterings out of the beam-stops, is observed for k < 5nm~^ indicative of an 
important degree of disorder. However, weak localized spots are apparent on the first 
one^l which are not visible on the second^. These spots reveals that an organization 
with some degree of lateral order on distances larger than the intermolecular one 
is possible, but their absence in the second case suggests that this organization can 
break down easily. These two extreme spectra most likely illustrate the sensitivity 
of the organization to the conditions of preparation and observation of the samples. 
They also show a subtle interplay between order and disorder whose approach can 
be described in two steps schematically. 

First, in order to obtain scattering in the wave vector range k < 5nm~^, peri- 
odic lattices with large parameters respecting a dense packing of the triple helices 
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1.3 nm 




Fig. 2. Left-handed collagen molecules (a) assemble in right-handed triple helices {b,c) repre- 
sented by rods (d) which in turn are proposed to be assembled with a regular shift (e) in order to 
create zones of gap " g " and overlap " o " giving account of the striations. 



close to that of a hexagonal lattice were proposed. For instance periodic defor- 
mations of a hexagonal lattice into a "quasi hexagonal" lattice with monoclinic 
symmetrj0 or a square-triangle lattice^ inspired by the crystalhne structures of 
model polypeptideJ^ were proposed. Second, in order to give account of the broad- 
ening of the discernable spots and of the presence of a background, some type of 
disorder was introduced in the organization. This was done considering a quasicrys- 
talline square-triangle tiling-'^- or paracrystalline fluctuations in a hexagonal lattice 
when the diffuse background is dominant!^. 

It must be noted that in all cases except the last the presence of a lattice should 
manifest itself as facets on the surface of the fibril^ which are never observed 



''This square-triangle lattice is usually labeled (3'^434) a notation pointing out the number of 
polygons around each of its vertices. 

'^Such facets are quite visible on toroidal aggregates resulting from DNA condensation, aggregates 




Fig. 3. Typical X ray pattern from Hulmes (a) accompanied by a drawing (b) strongly empha- 
sizing the observed scattering (for instance see also Wess and Orgel) and another pattern from 
Doucet obtained under different conditions of preparation of the sample (c) . The scattering along 
the equatorial direction is observed for wave vector A; < 5 nm~^ only, a fan-like background is 
always present. The scattering along the meridional direction is due to the longitudinal striations. 
Here k is defined by fc = 47rsinS/A. 



whatever the diameter of the fibriliP. The ir circ ular section was indeed emphasized 
as a major morphological fact quite earljUJJl^. propose here to consider that 
triple helices in fibrils adopt a lateral organization described by the algorithm of 
phyllotaxis, well known from botanists as a self-o rganized g rowth process and whose 
mathematical background is now well establishe d I I I This algorithm might be 



in which DNA is organized according to a hexagonal latticJEE^HSl. 
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well adapted here insofar as it leads to organizations with cylindrical symmetry in 
which each element has the most homogeneous and isotropic environment, in the 
same spirit as in ref. 1191 After having recalled the basis of phyllotaxis and made use 
of its algorithm as a template for overlap and gap regions of a fibril, we examine 
how the scattered intensities obtained from them can contribute to the analysis of 
the scattering data and to what extent phyllotaxis could provide a basis for the 
development of a structural model. 

2. Phyllotaxis 

A phyllotactic organization of points indexed by s is described by an algorithm 
such that the position of point s is given by its polar coordinates r = a^/s and 
6 = 27rAs that is r = {a/ ^/2tt\)\/9 which is the equation of a Fermat spiral here 
called the generative spiral. As the circle of radius r and area Trr^ = iTa?s contains 
s points, the area per point has the value no? (indeed it oscillates close to this 
value for small s and then converges towards it). It has been showrPQl that the 
the most homogeneous and isotropic environment, or the best packing efhciency in 
cylindrical symmetry, can be obtained for A = l/cj), where <j) is the irrational golden 
raticH (1 + -\/5)/2. The phyllotactic organization of a set of 2500 points with their 
Voronoi cells is shown on figure 

The generative spiral is scarcely discernible on this figure, the radius vectors of 
two consecutive points being at an angle 27r0 modulo 27r or —137.5 degrees, but 
the process generates quite conspicuous spirals, the so-called parastichies. There 
are three families of such spirals, each going through two opposite sides of a hexag- 
onal Voronoi cell, whose representatives for two of them with opposite chiralities 
are drawn on the figure. The code chosen for the Voronoi cells exemplifies their 
remarkable degree of organization. Pentagons and heptagons are concentrated in 
narrow circular rings with constant width which separate large rings of hexagons 
whose width increases progressively as one moves from the core to the periphery of 
the pattern. In the narrow rings pentagons and heptagons are associated in dipoles 
separated by hexagons whose shape is close to that of a square which would have 
two corners cut. From ring to ring, these dipoles are alternatively aligned along 
the spirals of one family of parastichies or the other. The counts of the numbers 
of dipoles and hexagons in each narrow ring and of the numbers of parastichies in 
each large ring, collected in table 1, put forward the mathematical substratum of 
the organization. All the numbers listed in the table are indeed members of the 
Fibonacci series, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, . . ., for which each number of rank u 
is defined as fu ~ fu-i + fu~2 from /o = 0, fi — 1. The number of sites in a narrow 
ring being that of its polygons, this number is also a Fibonacci number so that the 
radii of the narrow rings and the widths of the large rings are determined by the 

"^The two values A = l/(/>orA = l/(/>^ give two mirror images of the same points distribution as 
angles have to be considered modulo 2n and because 0^ = cf> + 1. 
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Fig. 4. A set of 2500 points organized according to tlie algorithm of phyllotaxis with the golden 
ratio. Each point is surrounded by its polygonal Voronoi cell whose number of sides corresponds to 
the number of first neighbors around the point. Blue cells are pentagons, red hexagons and green 
cells are heptagons. Those polygons are not necessarily regular and pentagons and heptagons are 
topological defects in the assembly of hexagons. Inserted on the left is a part of the figure showing 
in detail the evolution of Voronoi cells when crossing a ring of such defects. Representatives of two 
families of parastichies spirals are drawn, (a) being right-handed and (b) left-handed by convention. 
Moving from the core towards the periphery, new parastichies arc introduced when crossing a circle 
of defects. 



Fibonacci series. Moreover the distribution of di poles on each ring is typical of that 
of approximants of unidimensional quasicrystalJ^lES] 

One observe also that when a narrow ring is crossed, moving from one large ring 
to the following, the number of parastichies increases for one family only by exactly 
the number of dipoles contained in the narrow ring. Thus, each dipole aligned along 
one family of parastichies generate one more parastichy in this family, an effect 
quite equivalent to that of a dislocation adding one more reticular plane in part 
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Table 1. When going down the table, one moves from 
the core to the periphery of figure |4] All numbers ap- 
pearing in this table are successive Fibonacci numbers. 
This table stays limited to the beginning of the series. 



Narrow rings 


Large rings 


Dipoles of 
5 and 7-gons 


hexagons 


Right handed 
parastichies 


Left handed 
parastichies 


13 


8 










13 


21 


21 


13 










34 


21 


55 


34 










34 


55 


55 


34 










89 


55 


89 


55 








Fig. 5. One pentagon and one heptagon associated in a dipole create a dislocation in a bi- 
dimensional crystal of hexagons. The hexagons being kept regular the crystal is deformed by a 
step. 



of a crystal and represented in figure [5] It was indeed proposed to extend here tlie 
language of metallurgy, speaking of the large rings as hexagonal grains and of the 
narrow rings separating them as grain boundarieJ^H. The role of those dipoles is 
essential to ensure the best packing efficiency as the points of the new parastichies 
they introduce enable the Voronoi cells to maintain their area and keep the most 
isotropic shape possible. 

Thus, the local configuration is reproduced from one large ring to the next as 
it results from approxima nts in self-similar structures, this would not be the case 
if A was a rational 

numbeJSlEni. 

Such an organization must be understood consid- 
ering two kinds of disorders: a metric disorder corresponding to the fluctuations of 
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Fig. 6. Pattern with N = 2500 points obtained when one point over five along the generative 
spiral represented by small circles is suppressed in order to simulate vacancies in a gap region. 

distances between points with six first neighbors, their hexagonal cells being not 
regular everywhere, and a topological disorder associated to the introduction of 
pentagonal and heptagonal cells associated in dipoles. Those defects build concen- 
tric quasicrystalline rings with radii determined by the Fibonacci series so that an 
order is superposed onto the disordered lattice. This ordering within a disordered 
structure corresponds to the fact that, as the phyllotactic organization is defined 
by a single equation, it has, from the point of view of information theory, a low 
entropjO. 




3. A phyllotactic template for collagen fibrils 

We consider a fibril as an assembly of parallel rods, the triple helices, organized 
following the algorithm of phyllotaxis in the normal section of the fibril, but this 
model must include the fact that the triple helices do not build continuous straight 
rodH. As shown in figure[2t, the presence of striations along the fibrils led to propose 

"^It is remarkable that this solution appears in botanical organizations not because a long evolu- 
tionary process has optimized the structures but because it is the only one possible with respect 
to sequential growth in a self organizing process under geometrical constraints. 
'A more realistic representation of triple helices in fibrils should include the presence of an eventual 
torsion as suggested by several observations, however it is expected to be weak and should manifest 
itself out of the equatorial trace. 
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Fig. 7. Pattern with A'^ = 2500 points obtained when one point over five within close groups 
of five points, represented by small circles, is suppressed in order to simulate vacancies in a gap 
region. 



that the intermolecular interactions of the triple hehces put them in register so 
that each hne is interrupted on 35 nm every 300 nm, those interruptions being 
regularly shifted by 35 nm from line to line. These interruptions correspond to 
vacancies within the assembly of triple helices. A fibril can therefore be modeled 
as an alternate stacking of small dense cylinders containing triple helix segments, 
overlap regions, and less dense small cylinders containing four triple helix segments 
out of five only, gap regions, all with a thickness of 35 nm. 

For the cylinders corresponding to overlap regions, the distribution of the triple 
helices in their normal plane is that of the phyllotactic pattern of figure |4j without 
modification. For the cylinders corresponding to the gap regions, we assume a similar 
pattern, but with missing points or vacancies owing to the fact that one point over 
five must be discarded. We examine here two possible distributions of vacancies in a 
gap region which can be considered as two extreme cases as far as the homogeneity 
of the distribution is concerned. 

In the first case, the pattern is built eliminating one point over five along the 
generative spiral of the phyllotactic process. The result of this elimination is pre- 
sented in figure [6] in which eleven spirals of vacancies are separated by four spirals 
of points, a periodicity to be put in correspondence with the transverse repetition 
of figure Et. 

In the second case, the pattern is built eliminating one point over five within 
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Fig. 8. Radial average of the scattered intensity obtained from the phyllotactic pattern with 
A'^ = 6200 points, similar to that figure |4] overlap region, when its points are changed into disks 
having the diameter of the triple helices (1 nm), the form factor of the disk is the curve in grey. 

groups of five points determined by a pliyllotactic distribution of cells having an 
area five times larger than that of one point. This association of five triple helices 
is reminiscent of the microfibril model, but in which the microfibrils would have a 
phyllotactic organization. The result is presented in figure [7] in which the vacancies 
are distributed in a much more homogeneous manner than in figure IHl 

4. Intensity scattered by a phyllotactic organization of triple 
helices 

As the triple helices are perpendicular to the phyllotactic patterns, modified or not 
by the presence of vacancies, the equatorial traces of the intensities scattered by 
the cylinders, gaps or overlaps, containing them are those of these patterns. The 
scattered intensity expected from a distribution of points is the square of the am- 
plitude of its Fourier transform, but, as the points are representative of rods having 
a diameter of 1 nm, the points are to be changed into disks and the intensity de- 
creases as the form factor of the disk. The results are shown on figure 8, for the 
overlap region without vacancy, figure 9, for a gap region whose vacancies are orga- 
nized along spirals, and figure 10, for a gap region whose vacancies are distributed 
homogeneously. 

As for the spectra of figure |3l the scattering is concentrated in the wave vector 
range k < 5 nm^^, the absence of scattering for higher wave vectors is obviously due 
to the decrease of the form factor of the rod. Following an argument already used 
above, the quite intense and rather broad peak obtained for fc w 5 nm~^ is indicative 
of a disordered local packing. For fc < 5 nm^^, weaker structures are visible on the 
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Fig. 9. Radial average of the scattered intensity obtained from the pattern similar of that of 
figure|5](gap region with spiraling vacancies), with 6200 points instead of 2500. The arbitrary scale 
for intensities has been chosen to amplify small features at small k. 




Fig. 10. Radial average of the scattered intensity obtained from a pattern similar to that of 
figure [7] (gap region with homogeneously distributed vacancies), with 6200 points instead of 2500. 

three traces which are to be related to the particular mode of organization built 
by the algorithm of phyllotaxis, i.e. the presence of rings of dipoles in the overlap 
as well as in the gap regions and that of vacancies in the gap region only. The 
first contribution can be identified considering figure 11 which shows the intensity 
scattered by a pattern in which all the points at the centers of hexagonal Voronoi 
cells outside the rings of dipoles are blanked out. 

One can infer from these figures that the total scattered intensity is the sum of 
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Fig. 11. Radial average of the intensity scattered by circles of dipoles only. The contribution of 
the long range correlations between these circles is quite visible in the range of lowest wave vectors 
only. 

several contributions. One around A: w 5 nm^^ which corresponds to the local or- 
ganization within the rings. A series of peaks for k < 2.7 nm^^ which corresponds 
to the long range correlations between these ringf|sf. Broader and more complex 
scatterings for 2.3 < fc < 5 nm~^ in figures 9 and 10 therefore corresponds to the 
organizations of vacancies in the gap region. The contributions to the scattered 
intensity expected from the phyllotactic model of a fibril can therefore be summa- 
rized: 

- the presence of rings of dipoles in the overlap and gap regions for k < 2.7 nm~^, 
where their correlations lead to rather narrow scatterings whose separation decreases 
as fc 0, 

- the presence of vacancies in the gap regions for 2.3 < fc < 5 nm~^, where their 
correlation lead to broader and more complex scatterings, 

- the intermolecular distance for fc sa 5 nm~^, where the width of the peak reveals 
its local fiuctuations. 

5. Conclusions 

In a scattering pattern, the details which can be observed at a given value of the 
wave vector k have their origin in the density fiuctuations within a mobile test 
box of size proportional to fc^^. As, due to the form factor of the triple helix, the 
observations are limited to the wave vector range < fc < 5 nm~^, the intensities 
observed for fc < 5 nm^^ corresponds to density fluctuations on a scale greatly 

^The scattering corresponding to the overall circular shape of the pattern should manifest itself, 
but at very small wave vectors, owing to the large value of its radius. 
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larger than the intermolecular distance and the phyllotactic model identify them 
as topological defects and vacancies. The data presented in figure [3)d could then be 
read in the following way: 

- the beamstop and the equatorial trace issued from it at small k mask or contain 
most of the contributions coming from the rings of dipoles and the two first spots 
observed can be put in correspondence, owing to their positions and narrow widths, 
with the two of these contributions obtained at larger k, 

- the third group of broader spots can be put in correspondence with the set of 
lines identified as coming from vacancies in the gap region, but with a preference 
for their spiral organization, 

- finally the last group of still broader and more intense spots reveals, as already 
said, the disorder of the local molecular packing. 

Those suggested correspondences between observed and calculated intensity 
curves should be more deeply investigated: 

- carrying on the development of the phyllotactic model considering noble numbers 

different from the golden ratio and the possibility of the melting of the distribution 
91 . 

of dipoles , testing other distributions of vacancies in the gap region, taking into 
account the presence of a torsion and its influence on the off-equatorial scattering, 

- clarifying the influence of the conditions of preparation and observation on the 
structure, analyzing the central scattering which might contain information about 
eventual rings of dipoles and, finally, studying different types of collagen, particu- 
larly those without striations which might bring new information about the presence 
or distribution of vacancies. 

Up to now, the only thing which can be said is that the phyllotactic approach 
proposes a reasonable structural model somewhere at midway between the two 
extreme cases represented by figures 3b and 3c. This approach makes use of a very 
simple algorithm, which might be thought too simple to give account of such a 
complex structural problem as far as the relation between order and disorder is 
concerned. However, the interest of this algorithm is that, in the template so built, 
a topological order and a molecular disorder coexist intrinsically. This template 
belongs to the category of some models quoted above which discard true crystalline 
structures. It is also supported by the fact that very different systems facing similar 
constraints, compactness and radial symmetry, adopt this solution spontaneously, 
as discussed below. 

6. Comment on the growth process 

The fact that the simple algorithm of phyllotaxis was mostly used to describe vegetal 
growth calls into question its implementation in the case of collagen as the growth 
process of the fibrils, which proceeds by addition of triple helices at their periphery, 
differs from the one at work in plants. 

When this algorithm gives account of spiral organizations in flowers, the growth 
results from the addition of successive florets at the center of the flowers, the oldest 
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are pushed aside by the new ones and organize themselves on the generative spiral at 
azimuthal angles of 137.5 degrees from each other. This was illustrated studying the 
distribution of ferrofluid droplets falling down in a silicone oil in presence of a vertical 
inhomogeneous magnetic fields with cylindrical symmetry. The droplets polarized by 
the field push aside those already fallen in the oil and the positions of the successive 
droplets in the liquid evolve collectively towards that of the phyllotactic distribution 
characterized by the golden ratio each time a new one falla^. Recently a similar 
behavior has been observed for bubbles on a circular water surfacj24!. 

The phyllotactic organization of flowers is certainly the most notorious refer- 
ence, but it was also proposed to call for such an organization in a quite different 
problem, that of the the formation of Benard-Marangoni convection cells in a liquid 
layer contained in a cylindrical container and submitted to a vertical temperature 
gradient. As the gradient increases a transition takes place at which cells of uniform 
size determined by hydrodynamics appear in the liquid layer. Pentagonal and hep- 
tagonal cells, most often associated in dipoles, are observed among the hexagonal 
ones commonly formed in non cylindrical geometry. Their appearance was analyzed 
assuming the building of a phyllotactic pattern in order to respect the requirements 
of homogeneity and isotropy of the cell s in cy lindrical symmetry, whatever their 
distance from the center of the containeJ^lEir 

These examples stress upon the importance of cylindrical symmetry and of the 
necessity for the elements of the organization to move with respect to each other. 
We therefore propose here that the fibrils have a certain internal fluidity, for their 
molecules to reorganize each time a new one aggregates, and that their surface 
tension imposes them a circular section. Fibrils should then be considered within 
the frame of soft condensed matter studies rather than that of molecular crystals 
studies. If this were so, whereas in the cases of flowers or droplets the growth is 
generated from the center and in that of the convection cells the organization builds 
itself at a collective transition, this growth of fibrils by addition at their surface 
could also provide a third example of the power of the algorithm of phyllotaxis in 
structural studies. 
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